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Bootstrap approximation for the exchange-correlation kernel of time-dependent 

density functional theory 
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A new parameter-free approximation for the exchange-correlation kernel /xc of time-dependent 
density functional theory is proposed. This kernel is expressed as an algorithm in which the exact 
Dyson equation for the response as well as a further approximate condition are solved together self- 
consistently leading to a simple parameter-free kernel. We apply this to the calculation of optical 
spectra for various small bandgap (Ge, Si, GaAs, AIN, Ti02, SiC), large bandgap (C, LiF, Ar, Ne) 
and magnetic (NiO) insulators. The calculated spectra are in very good agreement with experiment 
for this diverse set of materials, highlighting the universal applicability of the new kernel. 

PACS numbers: 



The ab-initio calculation of optical absorption spectra 
of nano-structures and solids is a formidable task. The 
current state-of-the-art is based on many-body pertur- 
bation theory. A typical calculation involves two distinct 
steps: First, the quasi-particle spectral density function 
is calculated using the GW approximation, yielding ac- 
curate electron removal and addition energies, and there- 
fore a good prediction for the fundamental gap[l . In the 
second step, the Bethe-Salpeter equation (BSE) is solved 
using the one-body Green's function obtained in the GW 
step. Resonances, corresponding to bound electron-hole 
pairs called excitons, vifhich have energies inside the gap, 
can then appear in the spectrum. The two step procedure 
described above is a well-established method for yielding 
macroscopic dielectric tensors which are generally in good 
agreement with experiment [^-Jj. Unfortunately, solving 
the BSE involves diagonalizing a large matrix which cou- 
ples different Bloch state fc-points. As a consequence, the 
method is computationally expensive. 

Time-dependent density functional theory 
(TDDFT) H], which extends density functional the- 
ory into the time domain, is another method able, 
in principle, to determine neutral excitations of 
a system. Although formally exact, the predic- 
tions of TDDFT are only as good as the approx- 
imation of the exchange-correlation (xc) kernel: 
/xc(r,r',i- t') = Sv^c{r,t)/Sp{r',t'), where v^c is the 
TD exchange-correlation potential and p is the TD 
density. There are several such approximate kernels in 
existence, the earliest of which is the adiabatic local 
density approximation (ALDA)[9 , where Wxc(r,i) is 
determined from the usual ground-state local density 
approximation (LDA), calculated instantaneously for 
p(r,t). In practice however, the macroscopic dielectric 
function calculated using this kernel has two well-known 
deficiencies: the quasi-particle gap is too small and the 
physics of the bound electron-hole pair is totally missing 
- in fact ALDA does not improve on the results obtained 
within the random phase approximation (RPA) which 
corresponds to the trivial kernel /xc = 0[Tn]. In the 



present work we concentrate on the second of these 
problems, namely the missing excitonic peak in the 
spectrum. There have been previous attempts to solve 
this problem[ll). and there exist kernels which correctly 
reproduce the peaks in the optical spectrum associated 
with bound excitons. The nano-quanta kernel by Sottile 
et aL[12], derived from the four-point Bethe-Salpeter 
kernel, is very accurate but has the drawback of being 
nearly as computationally demanding as solving the BSE 
itself. The long-range correction (LRC) kcrncl^l^l E] 
has a particularly simple form in reciprocal space, 
/xc = —a/q'^, which limits its conrputational cost. This 
kernel produces the desired excitonic peak, but depends 
on the choice of the parameter a, which turns out to 
be strongly material-dependent, thereby limiting the 
predictiveness of this approximation. In the present 
work we propose a new parameter-free approximation 
for /xc, and demonstrate that this kernel is as accurate 
as BSE with a computational cost of ALDA. 

The exact relationship between the dielectric function 
£ and the kernel /xc for a periodic solid can be written 
as 



^'(q)xo(q,^) 

1 - [^(q) + /xc(q, w)] xo(q, w) 



= 1 



(1) 



where v is the bare Coulomb potential, x is the full re- 
sponse function, and xo is the response function of the 
non-interacting Kohn-Sham system. All these quantities 
are matrices in the basis of reciprocal lattice vectors G. 
The bootstrap kernel is a frequency-independent approx- 
imation given by: 



/xT(q,^) 



£ ^(q,^ ^ o)-t'(q) 

£o(q,w = 0) - 1 



(2) 



where £o(q, cj) = 1 — w(q)xo(q, i^)- e^"^(q, w = 0) is de- 
termined self-consistently with Eq. (IT]). We note that 
although Eq. (IT]) is exact, it is useful only when either 
/xc or £ is given; if neither are available then obviously it 



cannot be used as a generating equation for both quan- 
tities. With the addition of the approximation given by 
Eq. ([2]) however, both /xc and e can be determined from 
knowledge of xo exclusively. The modus operandi for do- 
ing so, is to start by setting /xc — and then solving 
Eq. ^ to obtain e"^ This is then used in Eq. ^ 
to find a new /xc, and the procedure repeated until self- 
consistency between the two equations is achieved. The 
advantages of this form for the kernel is that 1. it has the 
correct 1/q^ behavior [21 ITS]: 2. as e improves from cycle 
to cycle so does /xc; 3. the computation cost is minimal, 
as the most expensive part is the calculation of xo, which 
needs to calculated only once; and 4. most importantly, 
no system-dependent external parameter is required. 

The Xo in Eqs. ([l]) and ^ are, in practice, calculated 
using an approximate ground state xc functional, such 
as the LDA. To overcome the shortcomings of such an 
approximation, we further replace the xo by a model re- 
sponse function Xm coming cither from scissors-corrected 
LDA, or from GW or from LDA-|-t/. This has the advan- 
tage that Xm, and consequently x, has the correct gap to 
begin with. From the formal point of view, this replace- 
ment amounts to approximating the TDDFT kernel by 



Xo(q,w) 



1 



+ /r(q) (3) 



Using the method outlined above, optical spectra for var- 
ious extended systems [IS] were calculated using the full- 
potential linearized augmented plane wave (FP-LAPW) 
method [T7], implemented within the Elk code [T^. Ex- 
cept for the case of solid Ar, a shifted fc-point mesh of 
15 X 15 X 15 is used to ensure convergence [19]. In the case 
of solid Ar a shifted mesh of 25 x 25 x 25 fc-points was 
required for convergence of the optical spectrum. All 
the calculations were performed by scissor shifting the 
ground-state Kohn-Sham eigenvalues to reproduce the 
experimental bandgap. 

Presented in Fig. fl] are the results for some small (Ge 
~ 0.67 eV) to medium (diamond '^ 5.47 eV) bandgap 
semiconductors. For comparison, experimental data as 
well as the RPA spectra are also plotted. The experi- 
mental data clearly show that all these materials have 
weakly bound excitons leading to a small shifting of the 
spectral weight to lower energies compared to RPA. The 
results from TDDFT with the new kernel exactly follow 
this trend and are in excellent overall agreement with 
experiment. 

For Ge the TDDFT results are only slightly different 
from the RPA values which themselves are in agreement 
with experiment. It is clear that for Ge the RPA is 
enough and /xc does not significantly improve on the re- 
sult. This is in complete contrast to the spectrum of 
Si where the spectral weight is redistributed and, cor- 
responding to experiment, the TDDFT results show an 
enhanced Ei peak. The height of the E2 peak remains 
marginally overestimated by TDDFT. This overestima- 
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FIG. 1: (Color online) Imaginary part of the dielectric tensor 
(£2) as function[l9| of photon energy (in eV). Experimental 
data are taken from the following sources: Ge from 1201 Si 
from EH and m GaAs froml^ diamond fromi^, AIN from 
EH and SiC from|26] 



tion is not particular to the present approximation for 
/xc, it is also a feature of the BSE-derived kernel |^. The 
dielectric function for GaAs is also in very good agree- 
ment with experiment - subtle features like the kink at 
4.25 eV is well captured by our bootstrap procedure. 

The second column of Fig. [l] contains results for 
medium bandgap insulators. In all these materials a sig- 
nificant redistribution of the spectral weight to lower pho- 
ton energies is observed. For diamond, the bootstrap pro- 
cedure correctly leads to an enhancement of the shoulder 
at low photon energies. The position of the main peak 
around 12 eV is shifted to lower energies, and the whole 
spectrum is in near perfect agreement with experiment. 
AIN is a paticularly interesting case, TDDFT shifts the 
spectral weight to lower energies, and although the height 
of the peak is too large, the agreement with experiment 
is considerably better than that obtained by the equiv- 
alent BSE calculation P7H^. For SiC the results show 
an improvement over the RPA spectrum, but the height 
of the main peak as well as the shoulder at 9 eV are over- 
estimated. This trend is also observed in previous BSE 
results [30] . 

A stringent test for any approximate xc-kernel is in 
its ability to treat materials with strongly bound exci- 
tons. In these cases a new resonant peak appears in 
the bandgap itself and represents the bound state of an 
electron-hole pair. Perhaps the most studied test case 
for this phenomenon is the ionic solid LiF. Other exci- 
tonic materials which have also attracted attention and 
are considered particularly difficult to treat are the noble 
gas solids. Plotted in the first column of Fig. [2] arc the 
results for three materials of this class: LiF, solid Ar and 
Ne. What is immediately clear is that the bootstrap pro- 
cedure, which gave only a slight shift of spectral weight 
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FIG. 2: (Color online) Imaginary part of the dielectric tensor 
(£2) as function[19J of photon energy (in eV). Experimental 
data are taken from the following sources: LiF from 1311 Ar 
and Ne from El NiO from 1551 and Ti02 from El In the inset 
a smaller broadening is used to resolve the peaks better [19] 



for Ge, now gives rise to an entirely new bound excitonic 
peak inside the gap in all three cases. The location of the 
peak, which corresponds to the excitonic binding energy, 
is also very well reproduced for all these materials. 

Despite a good overall agreement we find that for LiF 
the main peak at 12.5 eV is overestimated, and the peak 
at 14.3 eV appears as a hump in the TDDFT results. 
Nevertheless, it is encouraging to note that the BSE spec- 
trum, as well that obtained using the BSE-derived ker- 
nel |35j . include a spurious peak at around 21 eV which 
is absent in the present calculations. Noble gas solids 
have very weak band dispersion and polarazibility, which 
results in very strongly bound electron-hole pairs. In 
the case of solid Ar, one can observe a strongly localised 
Frenkel exciton[36] at about 12 eV and a Wannier exciton 
at about 14 eV. This physics is totally missing within the 
RPA. Remarkably though, the bootstrap procedure cap- 
tures both these excitons, although the Wannier exciton 
is suppressed (see inset). Exactly like in BSE and LRC 
calculations J3S], the Frenkel exciton is under-bound by 
0.7eV. Ne has a strongly bound Frenkel exciton and the 
present calculations capture the corresponding excitonic 
peak. Similar to the BSE results [37], the height of this 
peak is overestimated by our TDDFT calculations. 

The second column of Fig. [2] consists of some special 
cases - NiO has an anti-ferromagnetic ground state, and 
the LDA+C/ method is needed to obtain a physically rea- 
sonable band structure for this material. This material 
provides the bootstrap technique with a test of its va- 
lidity for magnetic materials and also with a check of its 
performance when the scissors-corrected LDA is replaced 
by LDA+U, where U is chosen to reproduce the exper- 
imental gap. It is clear from Fig. [2] that the bootstrap 
method leads once again to the correct excitonic bind- 



ing energy. The experimental data for NiO are rather 
old and substantially broadened [33], and, assuming the 
veracity of these data, both TDDFT and BSE[38] over- 
estimate the peak height. It is worth noting that the 
BSE spectrum is redshifted relative to experiment, while 
the TDDFT spectrum is correctly positioned. Results for 
the anatase phase of Ti02 are also presented in Fig. [5] 
This material is important for its industrial use in pho- 
tovoltaics and has been well characterized using the BSE 
and GW method [311 EH] as well as experiment '34^ . Ti02 
is a useful test for the bootstrap method due to its non- 
cubic unit cell, which leads to directional anisotropy in 
the optical spectrum. As can be seen in Fig. [2] the boot- 
strap method captures this anisotropy very well indeed. 
Even subtle features like the small shoulder at ^ 4 eV in 
the out-of-plane dielectric function, which is missing in 
the in-plane case, are well reproduced. We find that our 
peak heights are slightly overestimated, which is also the 
case for the BSE results l39l. 
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FIG. 3: (Color online) Real part of the dielectric tensor (ei) as 
function of photon energy (in eV) for Si, GaAs and diamond. 

It is also interesting to compare the real part of the 
dielectric function with available experimental data. Re- 
sults for Si, GaAs and diamond are presented in Fig. [3] 
In all three cases, TDDFT results are in excellent agree- 
ment with the experimental data. We note that in the 
low frequency regime (below 2eV) the TDDFT results 
for GaAs deviate from experiments in exactly the same 
manner as found for the LRC kernel 1301. 
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FIG. 4; (Color online) Imaginary part of the dielectric tensor 
(£2) as function of photon energy (in eV) for LiF and diamond. 
£2 is plotted at various steps of the self-consistency. 

The effect of the bootstrap procedure can be seen 
in Fig. |4] where £2 is potted at each step of the self- 
consistency. As specified, the starting point is /xc — 
which yields the RPA spectrum in the first iteration. 



From then on, £2 changes substantially until convergence, 
indicating the importance of the self-consistency. 

An interesting point to note about Eq. (pi) is that this 
form of /xc is related to the two existing TDDFT kernels 
which capture the excitonic physics. In particular it cor- 
responds to the LRC kernel when a = Ane^^ /{sq — 1), 
and e is determined self-consistently. Also, like for the 
LRC kernel, the G = G' = part of f^c is by far the 
most important contribution. This also explains the fast 
convergence of the optical spectra with respect to the 
reciprocal lattice vectors included in /xc. Comparing 
our kernel to the BSE-derived approximation, in both 
cases the /xc is proportional to the screened Coulomb 
matrix elements e~^v. Choosing /xc to be proportional 
to the screened Coulomb interaction was also exploited 
by Turkowski et al. (see Eq. 1 and 15 of Ref. |4T|) . 
There remain several interesting aspects of the boot- 
strap method to be explored in the future: 1. the fi- 
nite q version of this approximation to determine the 
energy loss spectrum; 2. performance of this approxima- 
tion for two-dimensional systems, like graphene sheets or 
even nanotubes where excitonic effects are particularly 
strong; 3. to study the possibility of making this func- 
tional frequency-dependent; and 4. to explore analogous 
construction of /xc for finite systems and study its per- 
formance for excitonic spectra in molecular aggregates. 

Thus we have demonstrated that the bootstrap proce- 
dure gives a parameter-free TDDFT kernel which yields 
very accurate optical spectra. The same functional which 
produces a small shift, relative to the RPA, in the absorp- 
tion edge of Ge, also generates an entirely new excitonic 
peak within the bandgap of LiF, Ar and Ne. This in- 
dicates that the bootstrap kernel has wide applicability. 
The same quality of results can also be obtained by solv- 
ing the BSE, but the present method is highly desirable 
from the point of view of computation effort. 



* Electronic address: 'sharrQa@nipi-halle. nipg.de| 
" L. Hedin, Phys. Rev. 139, 796 (1965). 

G. Onida, L. Reining, and A. Rubio, Rev. Mod. Phys. 

74, 601 (2002). 

G. Onida, L. Reining, R. W. Godby, R. D. Sole, and 

W. Andreoni, Phys. Rev. Lett. 75, 818 (1995). 

S. Albrecht, L. Reining, R. D. Sole, and G. Onida, Phys. 

Rev. B 80, 4510 (1998). 

L. X. Benedict, E. L. Shirley, and R. B. Bohn, Phys. Rev. 

B 57, R9385 (2000). 

M. Rohlfing and S. G. Louie, Phys. Rev. Lett. 81, 2312 

(1998). 

A. Marini, C. Hogan, m. Griining, and D. Varsano, 

Comp. Phys. Comm. 180, 1392 (2009). 

E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 

(1984). 

E. K. U. Gross, F. J. Dobson, and M. Petersilka, Density 

Functional Theory (Springer, New York) (and references 



[10 

[11 

[12 

[13 
[14 
[15 
[16 

[17: 

[18 
[19 

[20 
[21 

[22 
[23; 
[24 

[25 

[26 
[27 
[28 
[29 

[3o; 

[31 

[32 

[33 

[34 
[35 
[36 
[37 



therein) (1994). 

V. L Gavrilenko and F. Bechstedt, Phys. Rev. B 55, 4343 
(1997). 

S. Botti, A. Schindlmayr, R. D. Sole, and L. Rening, Rep. 
Prog. Phys. 70, 357 (2007). 

F. Sottile, V. Olevano, and L. Reining, Phys. Rev. Lett. 
91, 056402 (2003). 

L. Reining, V. Olevano, A. Rubio, and G. Onida, Phys. 
Rev. Lett. 88, 066404 (2002). 

S. Botti, A. Fourreau, F. Nguyen, Y. Renault, F. Sottile, 
and L. Reining, Phys. Rev. B 72, 125203 (2005). 
P. Ghosez, X. Gonze, and R. W. Godby, Phys. Rev. B 
56, 12811 (1997). 

following experimental lattice parameters (in a.u.) were 
used: Ge 10.69, Si 10.26, GaAs 10.6870, AIN 8.2760, 
SiC 8.2386, diamond 6.7274, LiF 7.60804, Ar 10.030, Ne 
8.440, NiO 7.893389 and Ti02 a=7.1511, c=8.9899. 
D. J. Singh, Planewaves Pseudopotentials and the LAPW 
Method, Klu wer Academ ic Publishers, Boston (1994). 
(2004), URL |http : //elk . sourcef orge . net 
The calcuated spectra was broadened by (in eV) : Ge 0.11, 
Si 0.05, GaAs 0.1, AIN 0.2, SiC 0.23, diamond 0.22, LiF 
0.16, Ar 0.2 (and 0.1 for inset), Ne 0.2, NiO 0.2 and Ti02 
0.16. 

D. E. Aspnes and A. A. Studna, Phys. Rev. B 27, 985 
(1983). 

P. Lautenschlager, M. Garriga, L. Vina, and M. Cardona, 
Phys. Rev. B 36, 4821 (1987). 

D. E. Aspnes and A. A. Studna, Phys. Rev. B 27, 985 
(1983). 

P. Lautenschlager, M. Garriga, S. Logothetidis, and 
M. Cardona, Phys. Rev. B 35, 9174 (1987). 
D. Edwards and H. Philipp, Handbook of Optical Con- 
stants of Solids, edited by E.D. Palik (Academic Press, 
Orlando, 1985). 

V. Cimalla, V. Lebedev, U. Kaiser, R. Goldhahn, C. Fo- 
erster, J. Pezoldt, and O. Ambacher, Phys. Status Solidi 
C 2, 2199 (2005). 

S. Logothetidis and J. Petalas, J. Appl. Phys. 80, 1768 
(1996). 

F. Bechstedt, K. Seino, P. H. Hahn, and W. G. Schmidt, 
Phys. Rev. B 72, 245114 (2005). 

P. H. Hahn, K. seino, W. G. Schmidt, J. Furthmiillcr, 
and F. Bechstedt, Phys. Stat. Sol. (b) 242, 2720 (2005). 
R. Laskowski and N. E. Christensen, Phys. Stat. Sol. (b) 
244, 17 (2007). 

S. Botti, F. Sottile, N. Vast, V. Olevano, L. Reining, 
H. Weissker, A. Rubio, G. Onida, R. D. Sole, and R. W. 
Godby, Phys. Rev. B 69, 155112 (2004). 
D. Roessler and W. Walker, J. Opt. Soc. Am. 57, 835 
(1967). 

B. Sonntag, in Rare Gas Solids, edited by M. L. Klein and 
J. A. Venables (Academic Press, London, 1976), Vol. II, 
pp. 10211117. 

R. J. Powell and W. E. Spicer, Phys. Rev. B 2, 2182 
(1970). 

N. Hosaka, T. Sekiya, C. Satoko, and S. Kurita, J. Phys. 
Soc. Jpn. 66, 877 (1997). 

A. Marini, R. D. Sole, and A. Rubio, Phys. Rev. Lett. 
91, 256402 (2003). 

F. Sottile, M. Marsili, V. Olevano, and L. Reining, Phys. 
Rev. B 76, R161103 (2007). 

S. Galamic-Mulaomerovic and C. H. Patterson, Phys. 
Rev. B 72, 035127 (2005). 



[38] C. Rodl, Elektronische und exzitonische Anregungen in (2010). 

magnetischen Isolatoren(page 74) (2010). [41] V. Turkowski, A. Leonardo, and C. A. Ullrich, Phys. Rev. 

[39] H. M. Lawler, J. J. Rehr, F. Vila, S. D. Dalosto, E. L. B 79, 233201 (2009). 

Shirley, and Z. H. Levin, Phys. Rev. B 78, 205108 (2008). 
[40] W. Kang and M. S. Hybertsen, Phys. Rev. B 82, 085203 



